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1. 

p Abstract. We present J=2-l, J=3-2, J=4-3 ^''CO and 492 GHz [CI] maps as well as J=2"l and J=3-2 

■ ^"^CO measurements of the central regions in the nearby Sc galaxies IC 342 and MafFei 2. In both galaxies, the 
' distribution of CO and [CI] is strongly concentrated towards the center. These centers harbour modest starbursts. 

CO , Both galaxies have nearly identical ^^CO transitional ratios but the relative intensities of their ^''CO and [CI] 
emission lines differ significantly. The observed sets of line intensities require modelling with a multi-component 

, ^ molecular gas. In both galaxies, a dense component must be present (n(H2) ~ 10* — lO^cm"^) with kinetic 

^ ! temperatures Tki„ = 10-20 K (IC 342) or 20-60 K (Maffei 2), as well as a less dense (IC 342: a few hundred 

' cm~'^ at most; Maffei 2: « 3 x 10'^ cm"'^) and hotter (Tkin ~ 100-150 K) component. In both galajcies, neutral 

^— ■») , and ionized atomic carbon amounts are between 1.5 and 2.5 times those of CO. In both starburst centers about 

. half to two thirds of the molecular gas mass is associated with the hot PDR phase. The center of IC 342 contains 

' within R — 0.25 kpc an (atomic and molecular) gas mass of 1 xlO^M© and a peak face-on gas mass density of 

fT^ , about 70 Mq pc~^. For Maffei 2 these numbers are less clearly defined, mainly because of uncertainties in its 

■ distance and carbon abundance. We find a gass mass A/gas > 0.5 x IO^Mq, and a peak face-on gas mass density 
of about 35 Mq pc~^. 
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, 1. Introduction invisible at optical wavelengths. Mainly for this reason, 

IL" ' the distance of MafFei 2 is still quite uncertain. However, 

_^ . Molecular gas is a major constituent of the interstellar ^^^^^y assumed that it exceeds that of IC 342. We 

k> , medmm m galaxies. Withm the mner kiloparsec, many ^^^^ ^^^^ assumed a distance of 2.7 Mpc which is in line 

J, . spiral galaxies exhibit a strong concentration of molecu- ^^^^ overall cluster value, but 50 per cent greater than 

. , lar gas towards their nucleus. We have conducted a pro- reasonably well estabHshed distance to IC 342. Both 

gramme to observe a number of nearby galaxies m various 342 and Maffei 2 have been studied well at longer 

CO transitions and in the 492 GHz ^F,-^Fo CI transition wavelengths, including those of the CO line. Inasmuch 

m order to determine the physical condition of such central 342 is, together with NGC 253 and M 82, one of 

molecular gas concentrations. Results on NGC 253 (Israel, strongest (sub)millimeter line emitters in the sky, it 

White & Baas 1995), NGC 7331 (Israel & Baas 1999), been observed frequently. In the J=l-0 transition of 

NGC 6946 and M 83 = NGC 5236 (Israel & Baas 2001) i^cq, it was one of the first galaxies mapped at « V reso- 

have aheady been published. In this paper, we present the ^^^-^^ ^^^^^-^ ^ Lo 1978; Rickard & Palmer 1981; Young 

results for the very nearby galaxies IC 342 and Maffei 2 ^ g^^^^j^ ^982). More extensive mapping of the ^^CO 

whose basic properties are summarized m Table 1. 13^0 isotopes in both the J=l-0 and J=2-l transi- 

Both are major members of the IC 342/Mafrei group, tions was presented by Eckart et al. (1990) and Xie et al. 

which is located at a distance of about 2 Mpc (Huchtmeier (1994) Emission in the J=3-2 transition was detected and 

et al, 2000 and references therein). This appears to be the mapped by Ho et al. (1987), Steppe et al. (1990) and Irwin 

galaxy group closest to the Local Group but it is, unfortu- & Avery (1992), while the ^^CO isotope was detected in 

nately, located in the sky very close to the Galactic plane, this transition by Wall & Jaffe (1990) and Mauersberger 

Consequently, its members suffer high foreground extinc- et al. (1999). Finally, the J=4-3 and, with a somewhat 

tion rendering several of them, including Maffei 2, all but uncertain calibration, the J = 6-5 ^^CO transitions were 

detected by Giisten et al. (1993) and Harris et al. (1991). 
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Table 2. Observations Log 
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Table 1. Galaxy parameters 





IC 342 


Maffei 2 


Typo" 


SABcd 


SABbc 


Optical Centre: 






R.A. (B1950)'' 
Decl.(B1950)'' 


03''41™58.6" 


02''38'"08.0'' 


-f-67°56'26" 


-F59°23'24" 


R.A. (J2000) 


03'' 46"* 49.7'' 


Q2h4im54gs 


Dccl.(J2000) 


-|-68°05'45" 


-h59°36'll" 


Radio Centre : 






R.A. (B1950)= 


03''41'"57.3" 


02''38'"08.4'' 


Decl.(B1950)= 


-h67°56'27" 


-h59°23'30" 


+35 kms-i 


-31 kms-i 


Inclination i"^ 


25° 


67° 


Position angle P** 


39° 


26° 


Distance D" 


1.8 Mpc 


2.7 Mpc 


Luminosity Lg 


2 xlO^" Lb0 


1 xlO^o Lb0 


Scale 


115 "/kpc 


77 "/kpc 



Notes to Tabic 1: 

" RSA (Sandage & Tammann 1987) 

Dressel & Condon (1976); 
= Hummel & Grave (1990); Turner & Ho (1994) 

Rots (1979); Newton (1980); Hurt et al. (1996) 
" McCall (1989) for IC 342; for Maffei 2 sec text 
^ Buta & McCall (1999), assuming Galactic foreground for 
IC 342 = 3.32 mag (Madore & Preedman (1992); for 
Maffei 2 an uncertain foreground As = 8 mag is assumed. 



The ^^00/ ^"^CO isotope ratio was determined to be > 30 
by Henkel et al. (1998). High-resolution aperture synthe- 
sis maps of IC 342 have been published in J=l-0 ^^CO 



and "CO by Ishizuki et al. (1990), Turner & Hurt (1992), 
Wright et al., (1993) and Sakamoto et al. (1999), and in 
J=2-l 12C0 and "CO by Turner et al., (1993) and Meier 
et al. (2000). 

Maffei 2 was first detected in -'^^CO by Rickard et 
al. (1977) and subsequently mapped at 23" resolution 
by Weliachew et al. (1988). Emission in the J=2-l and 
J=3-2 transitions was measured by Wall et al. (1993) and 
Mauersberger et al. (1999). Maffei 2 was mapped in J=l- 
i^co and "CO by Ishiguro et al. (1989) and Hurt & 
Turner (1991), in J=2-l ^^CO by Sargent et al. (1985) 
and in J=3-2 ^^co by Hurt et al. (1993). 

Although the molecular line emission from Maffei 2 is 
somewhat weaker than that from IC 342, both are fre- 
quently included in the same observing programs. Both 
have been observed in a great variety of other molecular 
species: CS (Mauersberger & Hcnkel, 1989; Mauersberger 
et al., 1989; Sage et al., 1990; Paglione et al., 1995), 
HCN and HCO+ (Nguyen-Q-Rieu et al., 1992; Downes 
et al., 1992; Jackson et al. 1995; Paglione et al. 1997), 
H"CO+ and N2H+ (Mauersberger & Henkel, 1991); HNC 
(Hiittcmcistcr et al., 1995), CN and HC3N (Henkel et 
al. (1988), HNCO (Nguyen-Q-Rieu et al. 1991), OCS 
(Mauersberger et al., 1995) as well as H2CO and CH3OH 
(Hiittemeistcr et al., 1997). In particular the discovery of 
strong NH3 emission, first from IC 342 (Ho et al. 1990 and 
references therein), then from Maffei 2 (Henkel et al. 2000; 
Takano et al. 2000) showed that a significant fraction of 
the molecular gas in the centers of these galaxies is both 
dense and hot. 
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Fig. 1. Full resolution emission-line spectra observed towards the centers of IC 342 and Maffei 2. Top row: IC 342; 

bottom row: Maffei 2. Columns from left to right: J=2-l CO, J=3 2 CO, J=4-3 CO, [CI]. Vertical scale is actually 
in Tmb- Whenever available, ^^CO profiles are shown as the lower of the two profiles in the appropriate ^^CO box. 
Their brightness temperatures have been multiplied by three, so that they are depicted on the same temperature scale 
as [CI] in the right hand boxes. 



2. Observations 

All observations described in this paper were carried out 
with the 15m James Clerk Maxwell Telescope (JCMT) on 
Maima Kca (Hawaii) ^. Details arc given in Tabic 2. Up 
to 1993, we used a 2048 channel AOS backend covering 
a band of 500 MHz (650 km s"^ at 230 GHz). After that 
year, the DAS digital autocorrclator system was used in 
bands of 500 and 750 MHz. Integration times given in 
Table 2 are typical values used in mapping; central posi- 
tions were usually observed more than once and thus gen- 
erally have significantly longer integration times. Values 
listed are on-|-off. When sufficient free baseline was avail- 
able, we subtracted second or third order baselines from 
the profiles. In all other cases, linear baseline corrections 
were applied. All spectra were scaled to a main-beam 
brightness temperature, Tmb = T^/rj^y^; relevant values 
for r]mb are given in Table 2. Spectra of the central posi- 
tions in both galaxies are shown in Fig. 1 and summarized 
in Table 3. In Tabic 2, we have also listed the parameters 
describing the various maps obtained. All maps are close 
to fully sampled. In all maps, the mapping grid was ro- 
tated by the angle given in Table 2 so that the Y axis coin- 
cided with the galaxy major axis. The velocity-integrated 
maps shown in Figs. 2 and 3 have been rotated back, so 
that north is (again) at top and the coordinates are right 

^ The James Clerk Maxwell Telescope is operated on a 
joint basis between the United Kingdom Particle Physics and 
Astrophysics Council (PPARC), the Netherlands Organisation 
for Scientific Research (NWO) and the National Research 
Council of Canada (NRC). 



ascension and declination. For IC 342, the map grid origin 
in the J=2-l and J=3-2 CO maps is within a few arcsec- 
onds from the optical centre listed in Table 1 , whereas the 
J=4-3 and [CI] maps find their origin close to the radio 
centre; the two sets of maps are thus offset from one an- 
other by Aa, AS = +9", -3". For Maffei 2, the grid origin 
is again close to the optical centre; the radio centre occurs 
in the maps nominally at Aa, AS = -0.8", +6". 

3. Results 

3.1. CO distribution 

In both galaxies, molecular gas as traced by CO is strongly 
concentrated towards the center. The central molecular 
source is contained with R < 250 pc in the case of IC 342, 
and R < 200 x D/2.7 pc in the case of Maffei 2, where D 
is the true distance of Maffei 2 in Mpc. 

CO line aperture synthesis maps of IC 342 show, at 
resolutions of 2"-4", considerable structure only hinted at 
in our maps (Fig. 2). Clearly, the central CO concentra- 
tion is not toroidal in shape. Rather, the maps published 
by Sakamoto et al. (1999) and Meier et al. (2000) show 
individual CO peaks (cloud complexes) distributed along 
a perspectively foreshortened double spiral - although the 
foreshortening seems more than expected for an inclina- 
tion of only 25°. The molecular gas enhancement of spiral 
arms ends at about R = 250 pc from the nucleus. The 
nucleus itself is located at a minimum in the CO distri- 
bution. Emission from the high-density tracer HCN of- 
fers a very similar appearance (Downes et al. 1992). Our 
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Fig. 2. Contour maps of emission from IC 342 integrated over the velocity range Vlsr = -100 to 200 kms~^. North 
is at top. Offsets arc marked in arcminutcs with respect to J20G0 coordinates cto = 03^46'"49/7, So = +68°05 45 . 
Left to right: CO J=2-l, CO J=3-2, CO J=4-3 and [CI]. Contour values are linear in / TmbdV. Contour steps are 
16 Kkms-i (2-1, 3-2, 4-3) and 5 Kkms'^ (CI) and start at step 1. 




Fig. 3. Contour maps of emission from Maffei 2 integrated over the velocity range Vlsr = -200 to +150 kms~^. North 
is at top. Offsets are marked in arcminutcs with respect to J2GGG coordinates Oo = 02'*41'"54.*6, (5o = +59°36 11 . 
Left to right: CO J— 2-1, CO J=3-2, CO J=4-3. Contour values are linear in J T^hdV. Contour steps are 25 Kkms"^ 
(2-1 and 3-2) and 35 Kkms-^ (4-3) and start at step 1. 



J=2-l ^^CO map shows an elongated source, in which 
a smoothed version of the high-resolution map by Meier 
et al. (2000) is easily recognized. Emission in all our CO 
maps, as well as the [CI] map peaks at the position coin- 
ciding with cloud complex B in the designation by Downes 
et al. (1992). This peak becomes more pronounced with 
increasing J level, and in [CI]. This is partly the conse- 
quence of increasing angular resolution, but we also note 
that the peak coincides with the strongest concentration 
of thermal radio continuum emission (Turner & Ho 1983), 
thus with the region containing the highest spatial density 
of early type stars. The peak line intensities in Table 3 and 
the line ratios in Table 4 refer to this maximum. We thus 
expect them to exhibit properties commensurate with a 
starburst environment, i.e. those of a photon-dominated 
region (PDR). 



The region corresponding to complexes C and D 
(Downes et al. 1992), referred to as 'Eastern Ridge' by 
Eckart et al. (1990), behaves differently. It is almost as 
pronounced as region B in the J=2-l transition, but fades 
rapidly at higher J levels. Its molecular gas must thus be 
cooler on the whole than that of region B. This is, in any 
case, consistent with the much less impressive appearance 
of this region in the radio continuum maps by Turner & 
Ho (1983). There appear to be relatively weak, secondary 
max;ima in the neutral carbon emission at the northern 
and southern edges of the main CO distribution, in addi- 
tion to the primary [CI] maximum at cloud B. 

The major-axis position- velocity map of IC 342 covers 
a relatively small velocity range, in accordance with the 
mostly face-on orientation of the galaxy. With increas- 
ing J level, hence increasing resolution and decreasing 
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Notes to Table 3: a. Beam centered on nucleus; b. Total central concentration 

Table 4. Integrated line ratios in the centres of IC 342 and Maffei 2 



Transitions 


IC 342 

Nucleus 


Total Center 


Maflfei 2 
Nucleus 


Total Center 


i^CO (l-0)/(2 1)" 
i^CO (3-2)/(2 1)'' 
i^CO (4-3)7(2-1)'' 


0.95±0.1 
0.73±0.10 
0.67±0.10 


0.9 
0.7 
0.3 


0.9±0.2 
0.75±0.13 
0.65±0.10 


0.8: 
0.7 
0.5 


i^CO/i^CO (1-0)^= 
i^CO/^CO (2-1)'' 
i^CO/ "CO (3-2)<* 


10.7±1.3 
7.2±1.5 
10.2±1.4 




8.6±1.1 
10.0±1.4 
12.4±1.7 




CI/CO(2-l)'' 
011/00(2-1)" 


0.16±0.03 
0.54 


0.06 


0.10±0.05 
0.31 





Notes: a. IC 342 ratio from Eckart et al. (1990), Xie et al. (1994) and Meier et al. (2000); Maffei 2 ratio from J=l-0 data by 
Weliachew et al. (1988); both ratios at 21" resolution, b. This Paper, JCMT at 21" resolution; c. Weighted mean based on ratios 
presented by Rickard & Blitz (1985); NRAO at 65" resolution; Young & Sanders (1986); FCRAO at 45" resolution; Weliachew 
et al. (1988); IRAM 30 m at 24" resolution; Sage & Isbell (1991); NRAO 12m at 57" resolution; Xie et al. (1994) FCRAO at 
45" resolution; Meier et al. (2000) OVRO at 4."5 resolution d. This Paper; JCMT at 14" resolution, e. From Crawford et al. 
(1985) and Stacey et al. (1991), KAO at 55" resolution; JCMT convolved to 55" resolution. 



beamsmearing, the rotation steepens. In the J=4-3 CO 
p-V map the central material covers a range of about 80 
km over about 20". The rotation in this region is char- 
acterized by a velocity gradient (cf. J=4-3 CO in Fig. 4 
and J=l CO in Fig. 4 of Sakamoto et al. 1999) dV/de = 
7.8 kms^^/" , corresponding to dV/dRsv 2.1 kms~^/pc 
in the plane of the galaxy. 

The structure of the central CO source in Maffei 2 
is more complex than Fig. 3 suggests. The galaxy is 
much more tilted than IC 342, and the integrated veloc- 
ity images, showing only a resolved but featureless cen- 
tral source, are misleading. Maps of peak temperature 
Tmh rather than integrated-velocity temperature maps 
show a distinct double-peaked source (see also Hurt et 
al. 1993). This is in sharp contrast to IC 342, where peak- 



temperature and integrated-temperature maps show only 
marginally different morphologies. The deceptively simple 
appearance of Maffei 2 in Fig. 3 is easily explained by a 
glance at the major axis position-velocity maps in Fig. 5. 
The single peak in the integrated map results from the 
line-of-sight superposition of two distinct velocity compo- 
nents. This is, in fact, the same situation that we found 
to apply to NGC 253 (Israel et al. 1995), and it is one 
to watch out for in all high-inclination galaxies. The CO 
source in Fig. 3 becomes more compact with increasing J 
level, i.e. increasing frequency. As can be seen in Fig. 5, 
this is mostly the consequence of increasing resolution and 
decreasing beamsmearing. In the peak-temperature maps 
(not shown) the separation of the two peaks is only 13" 
(corresponding to 170 pc). The CO distribution shown in 
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Fig. 4. Position- velocity maps of CO emission from IC 342 in position angle 20°. Left to right: CO J=2-l, CO J=3-2 
and CO J=4-3. Contour values are linear in Tmb- All maps are one beamwidth wide, have a velocity resolution of 7.5 
km and have contour steps of 5 K, starting at step 1. Vertical scale is Vlsr 




40 20 -20 -40 40 20 -20 -40 40 20 -20 -40 



Y offset(arcsec) from 59 23 24.00 Y offset(arcsec) from 59 23 23.98 Y offset(arcsec) from 59 23 24.00 

Fig. 5. Position-velocity maps of CO emission from Maffei 2 in position angle 25°. Left to right: CO J=2-l, CO 
,7=3 2 and CO J=4 3. Contour values arc linear in Tmb- All maps arc one beamwidth wide, have a velocity resolution 
of 10 km s^^ and have contour steps of 3 K, starting at step 1. Vertical scale is Vlsr 



Figs. 3 and 5 is characteristic for a toroidal distribution 
of molecular gas, with a 45 - 60 pc wide gap in the cen- 
ter. However, the 'S'-shapcd morphology of CO in high- 
resolution (6") aperture synthesis maps by Ishiguro et al. 
(1989) and Hurt & Turner (1991) strongly suggests greatly 
enhanced emission from molecular gas in inner spiral arms 
instead. 



The effects of bcamsmcaring are quite noticeable in 
the p-V maps in Fig. 5. The average velocity gradient 
increases from 4.5 kms^^/" in the J=2 1 CO map to 
about 18 kms'V" in the J=4-3 CO map. The latter 
value corresponds to dV/dR = 1.4 kms~^/pc, i.e. not 
much different from the gradient in IC 342 (or NGC 6946 
and M 83 for that matter - see Israel & Baas 2001). We 
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note that the CO distribution in Fig. 5 just fiUs the steep 
leading edge of the major axis HI distribution in Maffei 2 
(Fig. 9 in Hurt et al. 1996). 



3.2. Line ratios 



We have determined the intensity ratio of the observed 
transitions at the centers of both galaxies; all intensi- 
ties are normalized to those of the J— 2 1 ^^CO line. 
The columns in Table 4 marked 'Nucleus' refer to ratios 
in a 21" beam centered on the nucleus, with the excep- 
tion of the [Cn]/C0(2-1) ratio which was determined af- 
ter convolution of the J=2-l ^^CO measurements to the 
55" beam of the [CII] observations, and the J = 1 — 
"'^^CO/"'^^CO isotopic ratio which represents a weighted 
mean from a number of determinations at various resolu- 
tions. The [CII] intensities (Crawford et al. 1985; Stacey 
et al. 1991), were converted to velocity-integrated temper- 
atures to obtain the line ratios in Table 4. The columns 
marked 'Total Center' refer to the intensities of the cen- 
tral concentration integrated over the source extent in the 
maps. The apparent decrease in source extent going from 
the lower J levels (lower frequencies) to the higher fre- 
quencies is mostly an artifact of limited resolution. When 
corrected for finite beamwidth, source dimensions at e.g. 
the J=2-l and J=4-3 transitions are in fact very sim- 
ilar. Nevertheless, the smaller area coverage in the [CI] 
and J=4-3 CO maps may ahve led to an underestimate 
of the total intensity. The entries in Table 4 suggest that 
this may indeed be the case. The J=l-0/J=2-l ratios 
have relatively large uncertainties, because we could not 
use observations of our own, but instead have used esti- 
mated J=l— intensities reduced to a 21" beam from the 
references given in the table. These ratios are, in any case, 
close to unity. 

The ^^CO ratios are quite similar for both galaxies, 
but the intensity of the [CI] and [CII] lines relative to CO 
is significantly lower in Maffei 2 than in IC 342, suggest- 
ing stronger PDR effects in the latter. The intensity of the 
[CII] line indicates the presence of both high temperatures 
and high gas densities in the medium as the critical values 
for this transition are Tkin > 91 K and n > 3500 cm~^. 
At the same time, such values must be reconciled with 
the much lower temperatures and (column) densities im- 
plied by the more modest isotopic ratios in the lower CO 
transitions. These isotopic ratios differ for the two galax- 
ies. For each individual transition the isotopic ratios just 
agree within the errors, but the overall patterns are clearly 
different. For Maffei 2, the isotopic ratio increases with 
increasing J level, but for IC 342 the isotopic ratio is min- 
imal in the J=2-l transition. These differences suggest 
significant differences in the molecular gas properties at 
the positions sampled in both galaxies 



4. Analysis 

4.1. Modelling of CO 

The observed ^^CO and ^^CO line intensities and ra- 
tios have been modelled with the large-velocity gradi- 
ent (LVG) radiative transfer models described by Jansen 
(1995) and Jansen et al. (1994). These provide model line 
intensities as a function of three input parameters: gas ki- 
netic temperature Ik, molecular hydrogen density n{H2) 
and CO column density per unit velocity {N{CO)/dV). 
By comparing model line ratios to the observed ratios 
we have identified the physical parameters best describ- 
ing the actual conditions in the observed source. Beam- 
averaged properties arc determined by comparing ob- 
served and model intensities. In principle, with seven mea- 
sured line intensities, properties of a single gas component 
arc overdctermined as only five independent observables 
are required. We found that no single-component fit could 
be made to the data of either IC 342 or Maffei 2. 

However, good fits based on two gas components can 

be obtained. In order to reduce the number of free pa- 
rameters, we assume identical CO isotopical abundances 
for both gas components and assign the specific value 
[i2CO]/[i3CO] = 40 (Mauersberger & Henkel 1993; Henkel 
et al. 1998). We identified acceptable fits by searching 
a grid of model parameter combinations (10 K < Tk < 
250 K, lO^cm--'* < ^(Ha) < lO'"^ cm--\ 6 x lO^^ cm'^ < 
N{CO)/dV < 3 X 10^^ cm"2) for model line ratios match- 
ing the observed set, with the relative contribution of the 
two components as the single free parameter. Largely as a 
consequence of the non-negligible finite errors in the line 
ratios, solutions are not unique, but rather delineate a 
range of values in a particular region of parameter space. 
To a certain extent, variations in input parameters may 
compensate for one another, leading to identical line ratios 
for somewhat different combinations of input parameters. 
We have rejected all solutions in which the denser gas com- 
ponent is also hotter than the more tenuous component, as 
we consider this physically unlikely on the linear scales ob- 
served. From the remainder of solutions, we have selected 
characteristic examples and listed these in Table 5. 

The line ratios observed for IC 342 allow only a rather 
limited range of gas parameters. One of the two compo- 
nents has well-determined parameters: a column density 
N{CO)/dV = 6 -10 xlO^^cm-^, an H2 density of 3000 
cm""^) and a kinetic temperature of 100 - 150 K. The 
other gas component either is less dense and equally hot, 
or denser and much cooler with a temperature of the or- 
der of 10 - 20 K. For Maffei 2, the choice of models is less 
restricted. Possible solutions pair a low-density, hot com- 
ponent and rather large column densities with a much 
denser, low-column density component of uncertain tem- 
perature. Alternatively, both components have a reason- 
ably high density of 3000 cm~^, but represent either a 
low temperature/high column-density or a high tempera- 
ture/low column-density gas. 
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Table 5. Model parameters 



Model Component 1 

Kinetic Gas CO Column Kinetic 

Temperature Density Density Temperature 

Tk n(H2) N{CO)/dV Tk 

(K) (cm-3) (cm-2(kms-i)-^) (K) 



Component 2 
Gas CO Column 

Density Density 
n(H2) N{CO)/dV 
(cm-3 (cm-2(kms-i)-i) 



Relative 
J = 2 - l^^CO Emission 
Component 1:2 



IC 342 



1 


150 


1000 


3x10" 


150 


3000 


0.6x10" 


0.65 : 


: 0.35 


2 


100 


1000 


3x10" 


100 


3000 


IxlO" 


0.35 : 


: 0.65 


3 


10 


10000 


3x10" 


100 


3000 


IxlO" 


0.25 : 


: 0.75 


Maffei 2 


4 


150 


500 


lOxlO" 


150 


3000 


0.3x10" 


0.25 


: 0.75 


5 


150 


100 


10x10" 


60 


10000 


0.6x10" 


0.50 : 


: 0.50 


6 


100 


100 


10x10" 


20 


100000 


0.3x10" 


0.30 


: 0.70 


7 


10 


3000 


6x10" 


150 


3000 


0.6x10" 


0.35 : 


: 0.65 



Table 6. Beam-averaged results 



Model Beam- Averaged Column Densities Total Central Mass Face-on Mass Density Relative Mass 
N{CO) N{C) N{R2) M(H2) Mga^ ^(Ha) agas Components 1:2 

(10i*cm-2) (10^1 cm-2) (IO^Mq) (Mo/pc-^) 









IC 342; Nh/Nc ■ 


= 2200; iV(HI)'" = 


1.5 X 10^° cm"^ 








1 


0.8 


2.3 


3.3 


0.5 


0.6 


49 


67 


0.7 : 


: 0.3 


2 


0.7 


3.8 


4.8 


0.7 


0.9 


70 


96 


0.4 : 


: 0.6 


3 


0.7 


1.0 


4.9 


0.7 


1.0 


72 


99 


0.3 : 


: 0.7 








Maffei 2; Na/Nc 


= 2500; iV(HI)'' = 


: 1.1 X 10^^ cm"^ 








4 


1.7 


1.6 


3.5 


0.6 


0.9 


28 


43 


0.55 


: 0.45 


5 


3.6 


1.3 


5.5 


0.9 


1.4 


43 


64 


0.8 


: 0.2 


6 


2.6 


1.1 


4.1 


0.7 


1.1 


32 


49 


0.7 : 


: 0.3 


7 


2.0 


1.9 


4.4 


0.7 


1.1 


32 


53 


0.55 : 


: 0.45 



Notes: a. Newton (1980); b. Hurt et al. (1996) 



A further check on the validity of the model results is 
provided by the C-'^^O measurements obtained by Eckart ct 
al. (1990). This isotope was not included in our modelling 
requirements. For isotopic ratios between ^^CO/C^^O = 
250 (model 1) and i^CO/Ci^^O = 225 (model 3), the pa- 
rameters of our solutions reproduce the observed intensity 
ratios in the J=l-0 and J=2-l transitions within 10%. 

4.2. Beam-averaged molecular gas properties 

The relation between carbon monoxide, neutral and ion- 
ized carbon combines the observed C and C+ intensities 
with the chemical models by Van Dishoeck & Black (1988) 
which show a strong dependence of the 7V(C)/A^(C0) col- 
umn density ratio on total carbon and molecular hydro- 
gen column densities. In the analysis of [CI], we took the 
kinetic temperatures, H2 densities and filling factors re- 
sulting from the CO analysis, and then solved for column 
density A^(CI). In practice, the column density of one com- 



ponent is frequently well-determined, whereas that of the 
other is more or less degenerate. For this reason, we solved 
for identical velocity dispersions in the two gas compo- 
nents. This procedure we also followed for [CII]; in those 
cases where both model kinetic, tc^iiipc^rrtturc; mikI gas den- 
sity are very different for the two components, we have 
modelled a single [CII] component with the higher kinetic 
tcmperatiire and the higher gas density of the two com- 
ponents. This applies specifically to models 5 and 6 (see 
below). Finally, we have related total carbon (i.e. C -I- CO) 
column densities to molecular hydrogen column densities 
by using an estimated [C]/[H] gas-phase abundance ratio. 

IC 342 has a measured central abundance [0]/[II] = 
2.0 X 10"^ (Vila-Costas & Edmunds 1992, Garnett et al. 
1997). We are unaware of any abundance determination 
for Maffei 2, and therefore use the mean value for the spi- 
ral galaxies listed by Zaritzky et al. (1994), [0]/[H] = 1.7 
X 10~^. Then, using the results obtained by Garnett et al. 
(1999), notably their Figs. 4 and 6, we have estimated car- 
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boil abundances [C]/[H] = 1.7±0.5 x IQ-^ and 1.45±0.5 
X 10~^ for IC 342 and Maffei 2 respectively. As a signifi- 
cant fraction of carbon is tied up in dust particles and thus 
unavailable in the gas-phase, we have adopted a fractional 
correction factor 5c = 0.27 (see for instance van Dishoeck 
& Black 1988), so that = [2iV(H2) + A^(HI)] = A 
[iV(CO) + A^(C)], where A = 2200 and 2500 for IC 342 
and Maffei 2 respectively. The numerical parameter A is 
uncertain by about a factor of two. 

In Table 6, we present beam-averaged column densi- 
ties for both CO and C (C° and C+), and II2 column 
densities derived under the assumptions just discussed. 
We also present the total masses estimated to be present 
in the central molecular concentration calculated from 
Ltot listed in Table 3, and the face-on mass densities im- 
plied by hydrogen column density and the galaxy inclina- 
tion. Notwithstanding the differences between the model 
cloud parameters, the beam-averaged results in Table 6 
are rather similar: hydrogen column densities, masses and 
mass-densities are well within a factor of two from one an- 
other. Beam-averaged neutral carbon to carbon monoxide 
column density ratios range from N{C°)/ N{CO) = 0.5- 
0.6 for IC 342 to iV(C°)/iV(CO) = 0.15-0.3 for Maffei 2. 
These values are well within the range found for other 
late-type galaxies such ad M 82, NGC 253, M 83, and 
NGC 6946 (White et al. 1994; Israel et al. 1995; Stutzki 
et al. 1997; Petitpas & Wilson 1998; Israel & Baas, 2001). 

4.3. The center of IC 342 

Even before the model fitting described in the previous 
section, the center of IC 342 was identified as the abode 
of hot and dense molecular gas, first by Martin et al. 
(1982) and Martin & Ho (1985) who found the presence 
of gas with kinetic temperatures estimated to be 70 K or 
higher. Somewhat later, measurements of the high-density 
tracer HCN led Downes et al. (1992) to suspect the pres- 
ence of molecular gas at various densities in the range 
"-(H2) = 150 — 3 X 10^ cm^"' and kinetic temperatures in 
the range Tkin = 50 70 K. This gas consisted of a few 
dense cloud complexes and much less dense gas inbctwecn. 
Multitransition HCN measurements (Jackson et al. 1995; 
Paglione et al. 1997) and other molecules (Paglione et al. 
1995; Mauersberger et al 1995; Hiittemeister et al 1997; 
Henkel et al 2000) suggest kinetic temperatures in the 
range of 50 - 200 K, densities of about 10'' cm~'^ as well as 
the presence of additional gas components at lower densi- 
ties. Such temperatures and densities are not unexpected 
for the site of a significant burst of star formation (cf. 
Turner & Ho, 1983; 1994). The PDR model calculations 
by Kaufman et al. (1999), applied to the relative intensi- 
ties of CO, [CI], [CII] and the far- infrared continuum of 
the center of IC 342 likewise suggest a temperature Tkin ~ 
150 K, a density around 10^ cm~^ and an ambient radia- 
tion field strength log Go = 2.5 (Go expressed in units of 
the Habing Field, i.e. 1.6 x 10~^ erg s~^ cm~^. All three 
model solutions in Table 5 are consistent with these val- 
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ues, although the higher densities suggested by the other 
molecular line observations are only provided by model 3. 
More support for model 3 follows from the C'^^O aper- 
ture synthesis observations carried out by Meier & Turner 
(2001) who find that the molecular cloud complexes have 
densities of the order of 1000 -3000 cm""^ and, more im- 
portantly, kinetic temperatures of 10-40 K (see also Meier 
et al. 2000). We thus accept model 3 as the most proba- 
ble approximation of the physical conditions in the central 
molecular gas of IC 342. About a third of the total molec- 
ular mass resides in a dense, fairly cold component with an 
excitation temperature close to the kinetic temperature, 
and very high optical depths in all observed CO transi- 
tions; even the ^^CO transitions are marginally optically 
thick with r = 1 2. Most of the remaining mass is hot, 
with Tex = 50 K. The J=l-0 ^^CO transition is opti- 
cally thin, but the higher transitions are optically thick 
although optical depths are not extremely high. 



Obviously, actual conditions will be more varied and 
complex, but even then we expect the beam-averaged re- 
sults to be not very different from the ones obtained here. 
We find a total central gas mass of about 10^ M0 which 
agrees well with estimates derived from e.g. HCN obser- 
vations, but is much lower than the mass values quoted 
in previous CO studies. However, these were all based on 
assumed CO to H2 ratios, notably the H2 column den- 
sity to CO intensity ratio X. The H2 column densities 
and masses presented in this paper are, in contrast, based 
on total gas-phase carbon amounts and the actual car- 
bon abundance. The uncertainty in these results is essen- 
tially determined by the error in the product 6c x [C]/[H], 
which is hard to quantify, but which we estimate to be 
about a factor of two. A mass of 10^ Mq is not very high, 
yet it is about 10% of the dynamical mass with R = 130 
pc; this renders any value much higher rather unlikely. 
Model 3 implies for IC 342 a value X = 3 x 10^^ H2 
mol cm~^/ K kms"'^, i.e. almost an order of magnitude 
lower than commonly used X values for the Galactic Solar 
Neighbourhood. However, such a low value is not uncom- 
mon for a galactic center. We have found similar values for 
NGC 6946 and M 83 (Israel & Baas 2001), and such values 
are by now solidly established for the center of our own 
Milky Way (Sodroski et al. 1995; Oka et al 1998; Dahmen 
et al. 1998). As an important consequence of these low val- 
ues for X and beam-averaged 7V(H2), the abundances of 
all molecular species that were derived using erroneously 
estimated much higher H2 column densities must be re- 
vised upwards by almost an order of magnitude. This is 
particularly important in the case of IC 342. Because it is 
one of the few galaxies bright in molecular line emission, 
it has extensively been used as the basis for abundance 
determinations of a relatively great variety of molecular 
species. 
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4.4. The center of MafFei 2 

In a similar vein, rough estimates of conditions to be 
found in the center of Maffei 2 fohow from observations of 
HCN, CS and other molecules: H2 densities in the range 
10^ — 10^ cm~'^ and probably closer to the former than to 
the latter, as well as kinetic temperatures of about 100 
K (Mauersbcrgcr ct al. 1989; Hiittemcistcr ct al. 1997; 
Paglione et al. 1997). The relative intensities of CO, CI, 
CII and the far-infrared continuum imply, within the con- 
text of the PDR model by Kaufman ct al. (1999), some- 
what higher densities close to lO^cm"^, temperatures of 
about 150 K and a radiation field log Go = 2.3. These 
temperatures and densities suggest that models 5 and 6 
are more appropriate than models 4 and 7 which lack the 
high densities that appear to be required. Models 5 and 6 
are not fundamentally different. They consist of a rather 
low density (n(H2) = 100 cm'^), hot (Tki„ = 100-150 K) 
and very optically thick (r = 15-20 for J=l-0) compo- 
nent and a more dense, cooler component of relatively low 
optical thickness. The actual temperature and density of 
the latter component is not firmly established. As can be 
seen in Table 5, temperatures between 20 and 60 K and 
densities between 10"^ and 10^ cm~^ can be obtained by 
varying the relative contributions of the two gas compo- 
nents, requiring only that the ratio T^[a/''T-i^2) is kept 
constant. The ionized carbon intensities can only be ex- 
plained by postulating the presence of a hot and dense gas 
phase not sampled by the CO observations but rather re- 
sulting from extensively photodissociated and ionized gas 
(interface with HII regions and supernova remnants?). A 
similar situation was found to apply in M 83 (Israel & 
Baas, 2001). 

In any case, the molecular medium in the center of 
Maffei 2 appears to be significantly different from that in 
the center of IC 342. In Maffei 2, the dense phase is warmer 
and probably denser, and its hot phase is much less dense. 
Maffei 2 is experiencing a strong central starburst compa- 
rable to the one in IC 342. A difference between the two, 
at least at radio wavelengths, is that Maffei 2 and M 83 
have much larger amounts of nonthermal emission than 
IC 342 (Turner & Ho 1983, 1994). This suggests that the 
starbursts in Maffei 2 and M 83 are more evolved than 
the one in IC 342, perhaps explaining the molecular gas 
differences as well. 

As the distance and metallicity of Maffei 2 have not 
been measured, we have had to assume plausible values. 
This is an additional source of uncertainty in the H2 col- 
umn densities, the masses and the mass densities derived 
for Maffei 2 in Table 6. As the observed HI column densi- 
ties are much less than the model-inferred total H column 
densities, H2 column densities, masses and mass densities 
scale linearly with the actual carbon abundance, whereas 
the H2 and gas masses scale with the actual distance 
squared. Because it is unlikely that the distance of Maffei 2 
differs much from the adopted value of 2.7 Mpc, the er- 
ror in the derived mass due to this uncertainty is not more 
than about a factor of two. The derived X value is subject 



to uncertainties in the carbon abundance but not to those 
in the distance. The N{ H2) values listed in Table 6 imply 
for Maffei 2 an H2/CO conversion ratio X = 2 — 3 x 10"^^ 
H2 mol cm~-^/ Kkrns"^, so that the same comment made 
for IC 342 also appUes to Maffei 2. 

5. Conclusions 

1. Maps of the central arcminute of the nearby starburst 
galaxies IC 342 and Maffei 2 in various transitions 
of ^^CO and ^^CO, and in [CI] confirm the compact 
nature of the central molecular gas emission in both 
galaxies. Most of this gas is within a 200 parsec from 
the nucleus. In both galaxies, the molecular gas seems 
to reside in bright inner spiral arms, rather than a disk 
or torus. 

2. Relative ^^CO line intensities observed in matched 
beams are virtually identical in IC 342 and Maffei 2. 
However, relative ^'^CO intensities, as well as the in- 
tensities of the [CI] and [CII] lines differ significantly, 
indicating different physical conditions for the molec- 
ular gas in the two galaxies. The observed line ratios 
require modelling with a multi-component molecular 
medium in both galaxies. 

3. In IC 342, a dense component with n{ H2) ~ 10^ cm~"^ 
and Tkin ~ 10 K is present together with a less dense 
n(H2) « 3 X lO^cm-^ and hotter Tk^n ~ 100 K com- 
ponent. Total carbon column densities are about 1.5 
times the CO column density. The modelling solution 
for IC 342 is reasonably well-defined. 

4. In Maffei 2, the parameters of the two components 
minimally required are less clearly defined. The more 
tenuous component has n{ H2) ~ 10^ cm~^ and Tkin ~ 
100 - 150 K, whereas the denser component has 
n( H2) « 10^ - 10^ cm-3 and Tkin = 20 - 60 K. Maffei 2 
seems to be more affected by CO dissociation as it has 
a C/CO ratio of two to three. 

5. In both starburst centers a significant fraction of the 
molecular mass (about half to two thirds) is associated 
with the hot PDR phase. 

6. With an estimated gas-phase [C]/[H] abundance of or- 
der 4 xlQ-'*, the centers of NGC 6946 and M 83 con- 
tain within R = 0.25 kpc similar (atomic and molec- 
ular) gas masses of about W'^Mq. Peak face-on gas 
mass densities arc typically 70 Mq pc^^ for IC 342 
and 35 M© pc~^ for Maffei 2, but the results for the 
latter are subject to relatively large errors caused by 
uncertainties in its distance and its carbon abundance. 
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